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1.  Introduction 


Consider  the  two-factor  nested  components-of-variance  model  given  by 
Yijk  "  w  +  Ai  +  BtJ  +  Cijk  for 

i  -  1,  2,  I  >  1;  j  -  1,  2,  ....  J  >  1;  and  k  -  1,  2 . K  >  1; 

where  E[A^]  =  0;  VarCA^]  =  02;  E[B^  ]  -  0;  Var[B^]  -  o|;  EtC^]  -  0; 

and  Var[C  ]  *  The  random  variables  Y  are  observable;  the 

ijk  c  ijk 

random  variables  A^,  . ..,  A^;  B^,  ...»  B^;  Cj^ . ^IJK  are 


pairwise  uncorrelated  and  unobservable  and  are  jointly  normally  distri¬ 
buted;  u,  02,  a2,  and  o 2  are  unobservable  parameters.  The  parameter 
space  Q  is  defined  by 


t(u» 


-  -  <  V  <  Cj2  *  o,  02  *  0,  02  *  0}. 


These  specifications  define  a  two-factor  nested  components-of-variance 
model  with  equal  numbers  in  the  subclasses  and  the  ANOVA  table  is  displayed 
in  Table  1. 

Table  1. 

ANOVA  table  for  two-factor  nested  components-of-variance  model 
with  equal  numbers  in  the  subclasses 


Source 

d.f. 

S.S. 

M.S. 

E.M.S. 

Total 

IJK 

£££I«k 

Mean 

1 

IJKY2 

Factor  A 

n^  -  1-1 

EES (Y.  -Y  )2 

1#  •  •  •  * 

Sl 

91  ■  °l «"!♦•*’* 

B  within  A 

n2  -  I(J-l) 

EEE  (Y.  ,  -T  )2 

1 J  •  1*  • 

S2 

«2  ’  “c^B 

Error 

n3  -  IJ(K-l) 

£££<Yuk-’ij.)J 

S3 

e3‘  °C 

In  this  model  there  are  several  functions  of  the  variance  components 


that  may  be  of  interest  in  applied  problems.  These  include 


°A*  °B*  0C*  0C/(°C  +  0B^  *  0C/(0C  +  °A)’  °A/  (°A  +  °B  +  °C) •  °B/' (°A  +  °B  +  °C)’ 

and  0^/(0 2  +  02  +  a2j#  The  oniy  functions  of  af,  a|,  given  above  for 
C  A  d  L  A  d  t* 

which  an  exact  size  confidence  interval  exists  is  and  0^/(0^  +  o|)  • 
Approximate  size  confidence  intervals  for  and  o|  have  been  given  by 
Moriguti  (1954),  Bulmer  (1956)  and  Howe  (1974).  Approximate  size  confidence 


intervals  for  o^J  (c*  +  +  a£)  t  a2/(02  +  o|  +  g2)  and  o*/ (p2  +  a2  +  02)  have 

been  given  by  Graybill  and  Wang  (1979).  In  this  paper  we  give  approximate  size 

confidence  intervals  for  o£/ (a^  +  a^) »  Q\!  +  0C^'  aA^aC  an<*  aA" 

Actually  we  obtain  approximate  size  confidence  intervals  for  o^/a^  only 

since  °a^°A  +  °C^’  an^  0C^0A  +  °C^  Can  Stained  from  these. 

In  Section  2  the  lower  limit  of  the  upper  confidence  interval  is  derived. 


in  Section  3  the  upper  limit  of  the  lower  confidence  interval  is  given,  and 

in  Section  4  is  a  short  discussion  of  other  methods  that  could  possibly  be 

used  for  confidence  intervals  on  o^/o^. 

A  C 


2.  Lower  Limit  of  the  Upper  Confidence  Interval  on  g2/p2 
Since  Y  ....  S2,  S2,  and  S2  are  complete  sufficient  statistics  for  this 

problem,  we  will  require  the  upper  confidence  interval  to  be  a  function  of  them. 
Write 

0,-0 

g(Y  ....  S2,  S2,  S2)  <  — —  <  - 

for  the  1  -  a  upper  confidence  interval  where  the  function  g(Y  s£,S2,S2), 

the  lower  confidence  point,  is  to  be  determined. 


6  -  6  JKa2 

Using  the  notation  in  Table  1  observe  that  - - -  =  — »  80  an  upper 


61  ‘  °2  63  °C 

confidence  interval  on  - - -  is  equivalent  to  an  upper  confidence  interval  on 


0. 


/Cf2. 


JA'~C 


Since  +  °q)  is  a  function  of  0^,  ©2,  03  only,  this  is  unchanged 


if  any  constant  c  is  added  to  Y^  ^  *n  the  model  given  in  Section  1.  Thus 


the  lower  confidence  point  g(Y  ...  S2,  S2,  S2)  should  also  be  unchanged  if  c 

thus  Y  ...  +  c  is  zero  and  S2,  S2,  S2 


is  added  to  Y  .  Let  c  *  -  Y  . 


are  unchanged  when  Y^^  I®  replaced  by  Yj^  +  c  (or  specifically  by  Y  ^  -  Y  ... 


Hence  g(Y  ...,  S2t  g2,  S2)  becomes  g(0,  S2,  S2,  S2)  and  the  lower  confidence 
point  is  a  function  of  S2,  S2,  and  S^  only.  So  the  objective  is  to  find  a 
function  of  S2,  S2,  S2,  say  f(S2,  S2,  S2)  such  that 

PCf(s2,  s2,  S2)  <  (01  -  02)/03] 

is  approximately  (and  very  close  to)  equal  to  a  specified  number  1  -  a. 


If  Y.^  is  rcplecc.1  by  cY  __  for  c  I1  0,  then  (6^  -  8 ^ ) /0 3  is 


unchanged.  Thus  we  require  f(c2S2,  czS2,  c2S2)  -  f(S2,  S2,  S2).  Let 

c2  -  1/S2,  then  f(S2,  S2,  s2)  -  f(S2/S2,  1,  S2/S2)  -  h(S2/S2,  S2/S2),  so 

the  lower  confidence  point  of  (0^  -  02>/03  is  a  function  of  S2/S^  and  S3/S^  only. 


-3- 


®1  ■  ®2 

So  a  large  sample  lower  confidence  point  for  - - -  is 


0  -  N  Xter(6) 
a 


s! 


0. 


'  NafclSl/S3  +  C2SS/S3  +  <C3S1/S3  '  %SJ/S3>  1 


ft 

S! 


g2  ^ 

—  -  1  -  N  I  {c. (S?/S?)  +  c,  +  (c.S2/S2  -  c,)2 
c2  0112  2  312  4 


S2 


-  q(S2/S2) 
?  12 


where  is  the  upper  a  probability  point  of  a  standard  normal  p.d.f. 

Therefore,  in  general  we  require  the  lower  confidence  point,  h(S2/S2,  S2/S2), 


61  "  02  S2 

of  — - -  to  be  0f  the  form  —  q(S2/S2),  and  we  determine  the  function 

^  n  2 


q(S2/S2)  such  that 


c2 

S3 


s2  22  91  -  ®2 
PC-2-  q(S2/S^)  <  -1- 

S2  1  2  e3 


(2.1) 


is  close  to  l-o.  We  require  q(S2/S2)  to  satisfy  (1),  (2),  (3)  below. 

(1)  When  the  hypothesis  H_:  a2  =  0  vs.  H  :  o2  >  0  is  accepted  for  a  size 

U  A  A  A 

a  test  the  confidence  Interval  should  include  zero,  and  when  is 

rejected,  h(S2/S2,  S2/S^)  should  be  an  increasing  function  of  S2/S2. 
To  test  HqJ  ■  0  vs.  Ha:  a2  >  0  the  hypothesis  Hq  is  accepted  if 

and  only  if  Sj/S2  <  Fa;  (This  test  is  uniformly  most  powerful 

unbiased).  Thus 

h(S2/S2,  S2/S?)  -  0  when  S  /S  <  F 

1232  12—  a* 

h(S^/S2,  S2/S^)  >  0  and  increasing  in  S2/S2  when  S2/S2  >  Fa;  n  q 
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g2 

Since  h(S^/S^,  S^/S^>  *  q(S^/s|)  we  obtain 

S3 


when  S^/S*  <  F 

1  2  -  o:  nltn2 


qCSj/Sp  -  0 

qCS^/S^)  >  0  and  Increasing  in  S^/S* 

(2)  When  J  -*•  “  (hence  -*■  »  and  n^  -*■  «)  the  confidence  interval  will  be 


when  SfVS*  >  F 

12  ot:  n^.n^ 


required  to  have  an  "exact"  confidence  coefficient  1  -  a.  When  J  ■+  00  it 
follows  that  00  and  n^  -*■  00  and  from  this  it  follows  that  S2  -*■  9^  in 

probability  and  S2  •+■  in  probability.  Start  with 


PC 


S2 

S1 


F 

a:n 


ll* 


<_  0^]  ■  1  -  a 


and  use  the  result  of  J  -*  «,  i.e.  replace  S^  and  S^  by  their  "equivalent" 


values  and  8^  respectively,  to  obtain 


S2  S2  9i  -  S', 

nJL  ( - i - i)  <  ■  ±  ■— Z]  -  1  -  a 


S2  S2F 
3  2  a:  n^, 


0. 


Hence  when  J  -»■  «> 


q(s2/s2)  -  0 


when  S?/S»  <  F 

1  2  —  a:  n^,00 


q(S?/S?)  -  - - - 1  when  S2/S2  >  F 

1  2  S2F  12  a: 


2  0:  n^,  ® 


0l-e2 


^  aA  <iu*ntity  — q -  is  dominated  by  and  we 

want 


A2  A2 

S2  S1 

P[—  x 


S2  S?F 

3  2  a:  n^»  n^ 


6ln 

- 

®3 


to  be  equal  to  1  -  a.  This  requires  q(S2/S2)  to  behave  like 

S2/S?F  for  large  S2/S2  in  the  sense  that 

i  l  a:  n^,  n^  12 


g* 

q(S?/s|)  -  - - -  (1  +  f.(S?/S?) }  where 


S?F 
2  a*,  n^,  n^ 


A(Sj/s|)  -*■  0  as  S2/S2  -v  » 

Any  function  q(s2/s2)  satisfying  conditions  (1),  (2),  (3)  will  give  an 
exact  confidence  coefficient  in  the  three  limiting  cases  0^/0 2  ■  1,  0^/02  *  " 
and  J  -*■<». . 

The  simplest  function  satisfying  those  conditions  is  the  linear  function 
qi(s2/s2)  ■  a^S^/S^  +  b^  where  a^  and  are  functions  of  n3,n2»n3,  and  a 

and  are  determined  by  the  conditions  (1),  (2),  and  (3).  However,  this  did  not 
give  results  as  good  as  desired  so  a  more  general  function  was  used,  namely 


q(S2/S2)  -  [aiS2/S2  +  bj_  +  C^Sj/S*)-1]/^.  n 


(2.2) 


From  condition  (3)  a  -  1. 

From  condition  (2)  b  (n, ,  «,«)--  F  i  c.  (n. ,  00 ,  00 )  *  0 

11  os  n^,  09  11 

From  condition  (1)  F  +  b-  +  c,/F  *  0  or 

a:  n^,  n^  1  1  a:  n^, 

c,  ■  -  F  (F  +  b-  ) . 

1  a:  n1,  n2'  a:  nj^,  n2  1 

Let  b^a^,  n2,  n^)  *  -  Fa.  n  w  for  all  n2  and  n3,  then 

c.  -  F  (F  -  F  ),  and 

1  a:  n^,  n2  a:  n^,  00  a:  n^,  n2 

q(S?/S?)  *  [S,2/S2  -  F  +  F  (F  -  F  )  S|/S?]/F  . 

mv  1  2'  1'  2  o:n  ,»  a:n^,n2  o:n^,»  a:n^,n2  2  1  arn^.n^ 

Thus  a  1  -  a  upper  confidence  interval  on  (9^-  02)/03  is  L2  <  (0^  “®2^®3  <  “ 

where  L2  is  defined  by 


l2  -  0 


if  S2/S2  <Fn,  n  „ 
12  a.  n^,  n^ 


(2.3) 


L*  ’  * T 2~<Sl'Sl  V  +  F«=nrn2  <F«=„1.-  -  Fa:Vn2>  *}'*¥ 


3  o:n^,n3 


if  S2/S2  >  F 


a:n3,n2 


-7- 


0 


Note  that  “  0  if  and  only  if  the  a  level  teat  of  H^:  o2  * 

is  accepted,  so  P[L.  =  0]  *  P[S2/S2  <  F  ]  <  1  -  a  and  P[L„  =  0]  =  1  -  a 

i  Li  —  a:  n^,  n2  —  *■ 

if  and  only  if  o2  «*  0.  The  probability  associated  with  Equation  (2.3)  is  a  function 
of  the  unknown  parameter  p  *  anc^  exact^y  equal  to  1  -  a  when  P  is  one 

or  infinity  or  when  J  is  infinite. 

The  excellence  of  this  approximation  is  indicated  by  Table  2, calculated  by 
simulation.  Columns  7,  8,  and  9  of  Table  2  contain  the  range  of  probabilities  of 
l<2  _<  (9^  -  as  t*ie  unknown  parameter  var*es  from  1  to  00 .  The 

approximation  appears  to  be  quite  satisfactory  even  for  small  sample  sizes. 

The  remainder  of  this  section  is  devoted  to  the  study  of  the  behavior  of 

52  0i  “  09 

P  >  p[-*-  q(S2/S2)  <  -=—r — -]  for  all  values  of  n. ,  n9  and  n-.  From  Table  2 

c2  1  i  -  0,  Li.  J 

53  J 

P  appears  to  get  closer  to  1  -  a  as  the  value  of  K  (hence  n^)  increases. 

In  fact  as  K  -*■  <*>  (hence  n^  -*•  °°)  the  problem  is  reduced  to  the  interval 

estimation  of  o2  in  the  one-factor  model  and  the  method  discussed  in  this 
A 

section  is  equivalent  to  Moriguti's  method  (1954).  From  this  one  knows  that  the 

„  _  2 

error  in  P  is  of  the  order  n22,  i.e.  P  ■  1  -  a  +  0(n2  )•  Another  way  to 

examine  the  behavior  of  P  is  to  expand  P  in  powers  of  n2  ^  and  n^  . 

The  algebraic  details  of  this  work  are  heavy  (see  Bulmer  (1957)).  The  resulting 
expansion  is 

P  -  1  -  o  +  <x0  +  a12/n2  +  a13/n3  +  a^/n2  +  a33/n3 
+  a23/n2n3  +  0(n“2,  n"3). 

This  assures  that  as  the  values  of  J  and  K  (hence  n2  and  n3)  increase 
the  accuracy  of  the  approximation  gets  better. 

In  Table  2,  P  is  between  0.9500  and  0.9597  when  I  -  3,  J  -  3,  K  -  3, 
and  1  -  a  -  0.95  and  when  I  -  7,  J  -  3,  K  -  3,  P  is  between  0.9500  and  0.9581. 
A  study  of  the  values  of  P  when  I  is  large  (n^  n2>  and  n3  are  large,  but 
^  -  n1/n2,  R2  -  n3/n3  remain  constant)  is  in  Wang  (1979). 
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3.  Upper  Limit  of  the  Lower  Confidence  Interval  on  ol/ai 


Since 


1  -  P[f(S2,S|,S2) 


we  use 


the  confidence  coefficient  a  in  the  lower  limit  of  the  upper  confidence 
interval  in  Equation  (2.3)  to  obtain  a  lower  1  -  a  confidence  interval 
on  (0^  -  02^®3  Siven  by  0<_  0^  -  d^)/  0^  <_  U  where 


U  = 


S3F1 


-ain^.n.^ 


CS1/S2  -  ri-iV. 


F  (¥ 

l-a:n^,n2  l-a:n 


1* 


-,l-..n1,«2)S2/Sl] 


if  S?/S?  >  F 

1  2  l-a:n  ,n2 


if  S?/s2  F 

1  2  l-a:n^,n2 


(3.1) 


We  could  determine  how  close  the  confidence  coefficient  of  this  confidence 
interval  is  to  the  nominal  1-a  by  simulation.  However,  due  to  the  expense 
of  computer  simulation  we  chose  a  different  route.  We  used  q^S^/S^)  =  aiSi^S2  +  bl 
and  conditions  similar  to  (1),  (2),  (3),  of  Section  2  to  obtain  the  confidence 
interval  where  is  given  by 


"i-0 


if  S?/s2  <  F, 


V  2  -  l-a:nlfn2 


(3.2) 


S2 


1  e2  2  1  -a:n1,n_  l-a:n.  ,n„^Fl-a:n  ,n„)  if  S^/S^  >  F 

S*  1’  3  12  1  3  12  1-a  :n,  ,n„ 

3  1  i 


Note  that  »  0  and  U  ■  0  if  and  only  if  the  1-a  level  test  of  =  0 

is  accepted.  Also  note  that  conditions  (2)  and  (3)  of  Section  2  are  satisfied 
by  the  confidence  Intervals  given  in  Equations  (3.1)  and  (3.2). 


i 


-9- 


--  -  'as-**** 


The  probability  associated  with  Equation  (3.2)  depends  on  the  value  of 


p  *  9i/02’  nl,n2’n3  anc*  can  ea8ily  calculated  if  n^  is  even;  we  get 
|  <S?/S2ri-a;n1.n3-ri-!i;„1,n2/Fl-<.:n1,„3>:l 


n„/2  ,  n_/2  nl/2_1 


-(TIT)  2  (d+T>  3  l  +  <ir>uv>y 


_  ,0y  c+1  2  vd+l'  3J 

y=0  y 1 2 


where  c  =  R. 


lFl-.:n1,n2/‘>>  d  ‘  <(,-1)RlR2Fl-<.:n1,»3/p  <s»e'*“>8  U979)). 


The  results  of  the  probabilities  of  (0 9 2 ) 3— ^1  are  ®^ven  ^able  3 

for  various  values  of  I,  J,  K  and  for  1-a  ■>  0.09,  0.95,  0.99.  The  actual 

probabilities  are  quite  close  to  the  specified  probabilities  even  for  small 
sample  sizes.  We  expect  the  results  to  be  even  better  if  the  more  general 
confidence  interval  0^(9^  “  02^®3  —  u*s  used  where  U  is  given  in  Equation 

(3.1). 


4.  Comparison  with  Other  Methods. 


The  literature  does  not  contain  any  references  that  have  been  evaluated  and 
directly  relate  to  confidence  intervals  on  o^/a£.  Perhaps  Satterthwaite's  (1946) 
method  could  be  used  but  this  procedure  is  extremely  poor  when  used  to  place 
confidence  intervals  on  the  difference  of  expected  mean  squares  (i.e.  on 
=  J  K  a^r/a?,  ).  Broemeling  (1969)  presents  a  method  for  placing  simultaneous 
confidence  intervals  on  and  0^/0^..  This  method  can  be  used  to  place  con- 

AO  D  0 

fidence  intervals  on  o^/o^. 


We  use  Equation  (15)  in  Broemeling  (1969)  to  obtain 
p[0<kjc|/o|  <  Sj/sj  r2_  ^  ^3  >  <1  -  <0* 


(4.1) 
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which  can  be  used  for  a  lower  confidence  interval  on  KJc2/o2  with  confidence 

A  C 

coefficient  greater  than  or  equal  to  (1  -  a)2.  Clearly  the  1  -  a  lower 
confidence  interval  in  Equation  (3.2)  above  is  shorter  than  the  (1  -  a)2 
confidence  interval  in  Equation  (4.1).  Thus  the  confidence  interval  on 
°A^°C  <*erived  from  the  procedure  by  Broemeling  is  not  as  good  as  the  method  pre¬ 
sented  in  this  paper. 
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Table 


1A  t-  t^CO  f- 
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